Phonon structure in I-V characteristic of MgB 2 point-contacts 
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The search of the phonon structure at the above-gap energies was carried out for d 2 V/dI 2 (V) 
spectra of MgB2 point contacts with a normal metal. The two-band model is assumed not only for 
the gap structure in dF/d/(F)-characteristics, but also for phonons in d 2 V/dI 2 (V) point-contact 
spectra, with up to the maximum lattice vibration energy. Since the current is carried mostly by 
charges of 3D-band, whereas the strong electron-phonon interaction occurs in 2D-band, we observe 
the phonon peculiarities due to "proximity" effect in fc-space, which depends both on the contact 
orientation with respect to the crystal axes and the variation of interband coupling through the 
elastic scattering. 

PACS numbers: 74.25Fy, 74.80.Fp, 73.40.Jn 



Introduction. 

The superconducting properties of the newly discov- 
ered s-p compound MgB2 have been attracting re- 
cently much attention of the superconductive commu- 
nity. The relatively high T c ~ 39 K, the absence of mag- 
netic order and strongly correlated charge carriers, on the 
one hand, and the strong links between grains for large 
enough critical currents and magnetic fields, on the other 
hand, make it very promising for complete understand- 
ing of the superconducting mechanism and its successful 
applications at moderate temperatures 

Unfortunately, the simplicity of this compound does 
not go too far for immediate understanding of all details 
of its superconducting mechanism. From the theoretical 
point of view, this compound appears to be a rare ex- 
ample of several (at least two) disconnected bands of the 
Fermi surface (FS) with quite different dynamical proper- 
ties ||, |[ |j| . The theory predicts that one of those is two- 
dimensional (2D), with strong electron-phonon interac- 
tion (EPI) , while the other is 3D with weak EPI and even 
larger density of states at the Fermi energy |l| ||, ^| . The 
interband scattering between these two parts is weak [ 5 1 , 
and they behave like an "intrinsic proximity effect" [7\. 
The superconducting and normal parts are not separated 
spatially, like in the McMillan model §, |] of SN sand- 
wiches, but are weakly connected in fc-space through the 
interband scattering. In connection with this behavior, 
one has to recall earlier works of two-band superconduc- 
tivity |h| O, 12 1 . But there is an important difference 



between the predictions of those old papers and what is 
observed in Ref. and some other point-contact and tun- 
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neling works (see, for example, Refs.g [L§). Since the 
2D and 3D parts of Fermi surface are very anisotropic, 
the tunneling and point-contact current depends strongly 
on the orientation with respect to the crystal axes, and 
for c-orientation the contribution of the 3D band prevails 
ifilH , Superconducting energy gap values, corresponding 
to these two bands, differ very much (up to the factor of 
3 Jl5|, [l6|]), and turn to zero at the same critical temper- 
ature ~ 38 — 40 K. Moreover, if the separate energy gaps 
converge to a single averaged gap, the critical tempera- 
ture remains almost the same. Possible explanation of 
this behavior |l5| is that the charges of the 3D band play 
much greater role in the PC current, due to the higher 
Fermi velocity and larger plasma frequency. By analogy, 
the situation looks as if one measures superconducting 
energy gap only from the "normal' side of the McMillan 
sandwiches in coordinate space. T c does not decrease no- 
ticeably with increasing elastic scattering, since the latter 
is not strong enough. 

Luckily, in MgB2 no other bosonic excitation, except 
phonons, is anticipated ||, [l(| Thus, the only mech- 
anism of the formation of Cooper pairs is EPI, although 
when a new superconductor with relatively high T c is 
discovered, other mechanisms are always proposed. 

In spite of great variety of spectroscopic measurements, 
which convincingly show the existence of two different su- 
perconducting energy gaps belonging to the above men- 
tioned two parts of FS (see, for example, Refs. 13, [l8| , 
and more recent papers [M H), there are only a few re- 
ports of direct experimental proof for phonon structure 
in the point-contact characteristic jl{J or tunneling |2t| . 

In this paper we present the first attempt to identify 
the structure of MgB2 point-contact (PC) spectra at en- 
ergies up to the maximum phonon frequency. We show 
that the observed peculiarities are unequivocally caused 
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by phonons, and their appearance correlates with the ob- 
served superconducting energy gap for a given point con- 
tact. In our opinion, this is due to the strong anisotropy 
of multi-band electron and phonon spectra, and to the 
influence of the elastic scattering, which determines both 
the intraband and interband EPI coupling in this com- 
pound. The random deviation of the contact axis from 
the predominantly c-direction plays a crucial role. 
Specific of experiment 

The samples are c-axis oriented films of MgB2, whose 
characteristics are described in Refs.|l8l Eli. The con- 



tacts are made by direct touch of MgB2 with a sharpened 
normal metal electrode (Cu or Ag) . The axis of the con- 
tact nominally coincides with the c-axis of the film, and 
thus, determines the main direction of the current. The 
orientation selectivity in the point-contact spectroscopy 
is low (lower, than in plane tunnel junction, where it goes 
exponentially), and is within a solid angle of roughly a 
few tens of degrees. Moreover, since not every touch pro- 
duces the desired appearance of the energy-gap I-V char- 
acteristic, we moved the electrodes relative to each other 
in order to penetrate through the accidentally damaged 
surface layer. This operation changes the contribution of 
the perpendicular aft-direction and produces additional 
defects in the contact region and, thus, shortens further 
the small electron mean free path. In a few experiments 
we tried to make a contact with a broken side-face of the 
film, in order to increase the occurrence of a6-direction 
for the current. 

Thus, our "input" characteristic is the shape of su- 
perconducting energy-gap structure, which is described 
elsewhere jl8|, The main goal is to investigate the 
first and second derivatives of the I-V characteristic at 
large biases, providing the contact survives the long-term 
measurements. 

The attractive feature of PC spectroscopy is that the 
measured second derivative for the normal state turns out 
to be directly proportional to EPI function a 2 (uj)F(ui) 
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where a, roughly speaking, describes the strength of the 
electron interaction with one or another phonon branch 
for electron transport through the contact with diame- 
ter d, F(ui) stands for the phonon density of states. In 
PC spectra the EPI spectral function is modified by the 
transport factor, which increases strongly the backscat- 
tering processes. 

However, in the superconducting state the above-gap 
singularity at phonon frequencies can be caused by sev- 
eral mechanisms 
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The features, which are 
strongly shifted along the voltage scale upon increas- 
ing the temperature and field, should be disregarded, as 
due to the nonequilibrium superconducting phenomena 



p4| , |25[ |. Only those, whose positions do not depend on 
T, H, except the small shifts of the order of supercon- 
ducting energy gaps (2-j-7 meV), are considered. 

The temperature is varied from 1.5 K up to T > T c 
with the magnetic field up to 4.5 T applied along the 
afe-plane. At low temperatures this field is not enough 
to destroy superconductivity. In order to observe the 
inelastic point-contact spectrum in the normal state, we 
raise both the temperature and the field. 

The magnetic field docs not influence the inelastic PC 
spectrum (1), because of the short mean free path. But 
it can strongly suppress superconductivity in the weakly 
superconducting 3D band, and destroy the exchange of 
the Cooper pairs from the strongly superconducting 2D- 
band. Since the current flows mostly from the 3D-band, 
the phonon structure appears indirectly, its intensity be- 
ing strongly influenced by the field. Only small admix- 
ture of the direct 2D-band contributes to EPI, due to 
random inclusion of the a&-direction in the current. 

Results and Discussion. 

In Fig. la the dV/dI(V) and V 2 [eV) cx d 2 V/dI 2 (eV) 
characteristics are shown in the superconducting state 
at zero field. In the energy-gap region the two gaps are 
clearly seen with Ai = 2.4 and A2 = 5.7 meV. Fur- 
ther on, we simply denote by A the position of dV/dl 
minimum. For many junctions we determine the su- 
perconducting energy gap by the BTK-fitting procedure 
|p"8f . For PC with not too large T-parameter, intro- 
duced by Dynes, this value does not differ too much 
from the true energy gap. The first derivative at above- 
gap energies is slightly asymmetric and nearly parabolic, 
which leads to an approximately linear background in 
d 2 V/dI 2 (V). This background is subtracted from the 
raw V2(V)-dependence, an d the result is shown by curve 
2 (panel a). The latter curve displays an antisymmetric 
structure, whose amplitude is about a fraction of percent 
in differential resistance. In the panel (b), this structure 
is compared with the phonon density of states (PDOS) 
p6[ . Good correspondence is seen between the positions 
of PDOS peaks and maxima in d 2 V/dI 2 (V). 

The reproducibility of this PC spectra for different con- 
tacts can be seen in Fig. 2a, b, where the raw data for 
three different junctions are shown. Here, the contact 
resistance varies from 43 to 111 51 with gap minima at 
A = 2.1, 2.6, and 4.86 meV, respectively. All the second 
derivative spectra (b) correlate with PDOS. The slight 
variation is probably due to the anisotropy and different 
scattering rates in the contacts. Compared with the the- 
oretical EPI function (see, for example, Ref. || 0), the 
peak at eV=60-^70 meV of Ei g boron mode is not too 
much higher in intensity, in accord with observation of 
Ref.@. 

According to our classifications of the energy gap struc- 
ture JlH , it is due to the random orientation of the con- 
tact axis and scattering of charge carriers between two 
band of Fermi surface, having gaps Ai and A2 mixed. 
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FIG. 1: a) Curve 1: Differential resistance R = dV/dI(V) 
for MgB2-Ag point contact. The structure at bias jlO mV 
shows two superconducting gaps corresponding to two differ- 
ent parts of the Fermi surface. The junction is made with 
the broken side-face of the film. Curve 2: The point-contact 
spectrum V2 oc d 2 V/dI 2 of the same junction (V2 is the second 
harmonic voltage with modulation voltage V\ = 2.52 mV). 
d{lnR)/dV = 2V2V2/V 2 , where R = dV/dI(V). The linear 
background is subtracted from the raw data. 
T = 4.2 K, H = 0, R = 35 . 

b) The point-contact spectrum averaged between plus and 
minus ^-polarity (the same, as in panel a), compared to the 
experimental phonon density of states of MgE>2 taken from 
Ref. ^6|. The ordinate axis for PDOS is given in arbitrary 
units. 



With an increase of interband scattering the magnitude 
of Ai (A2) moves to the higher (lower) value, respec- 
tively. For dirty contacts, where the admixture of the 
2D-band is essential, only one maximum remains with a 
broad distribution around ~ 3.5 meV Q, while T c re- 
mains almost the same. In this case I is smaller than 
d (where I and d are the electron mean free path and 
the size of the contact, respectively), and the inelastic 
backscattering contribution to the phonon structure, pro- 
portional to l/d p3|, is small. Therefore, in the super- 
conducting state, the observed phonon structure presents 
mainly the elastic contribution to the excess current [ p4| . 
The elastic term is proportional to the energy depen- 



dent part of the excess current I exc (eV), similar to the 
phonon structure in the quasiparticle DOS for tunneling 
spectroscopy p3, |2j| : 
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where A ira (eV) is the gap parameter in the 3D band, in- 
duced by 2D-band EPI. As seen in Fig. 2, if the A value 
increases (panel a, curves 1-3), then the intensity of the 
phonon structure in the units of d(ln R) / dV = 2\/2 V% /V^ 
increases too (panel b). That is because the modulation 
voltage V\ decreases, whereas the amplitudes in Vi[eV\ 
units remain approximately the same. Note, that for 
curve 1 in Fig. 2a one has to take the lower gap Ai, be- 
cause the current is mainly determined by 3D-band, as 
mentioned above. 

In Fig. 3, the normal state spectrum is displayed to- 
gether with the differential resistances in the supercon- 
ducting state, showing the energy gap structure. In spite 
of a 10 times increase of the temperature smearing (~ 20 
meV at 41 K instead of 2 meV at 4.2 K), the residual 
phonon structure is still visible in Fig. 3b. The smeared 
phonon features are superimposed on the rising linear 
background. For curve 2, we see an increase in scatter- 
ing at ~ 35 meV, where the acoustic phonon peak occurs, 
and the saturation at ~100 meV, where the phonon spec- 
trum ends. In the normal state, only those nonlinearities 



remain, which are due to the inelastic processes 22 . We 
stress that judging from the superconducting gap value 
(~ 3.69 meV for curve 2 in panel a) the essential contribu- 
tion is expected from the A 2 (E) of 2D-band for normal- 
state spectrum No. 2. This spectrum should contrast with 
curve 6 in Fig. 4 (see below), where the direct contribu- 
tion from the 2D-band is small, due to lower value of 
A. Thus, the backscattering processes from the 3D-band 
are mostly essential here. The same is true for curve 6 
in Fig. 4, where the phonon features are not seen at high 
energies. The shape of spectrum 1 in Fig. 3b presents the 
behavior intermediate between these two extremes, al- 
though its energy gap is approximately equal to curve 6 
in Fig. 4. Beside the well seen phonon spectral features in 
the range 30 -j- 100 meV, it possesses a low energy bump 
at about 20 meV. 

Note, that the modulation voltages for spectra 1,2 in 
Figs. 3b and 4 are the same. That means that they have 
the same second derivative scales, and the intensity of 
the second derivative is greater for curve 2 in Fig. 3b, 
than those of curve 1 in Fig. 3b and curve 6 in Fig. 4. 

The increasing second derivative of the TV character- 
istic in the normal state is due to excess generation of 
nonequilibrium phonons at the contact. The PC spectra 
with phonon features in the normal state at T> T c , like 
those shown in Fig. 3b, usually exhibit the nonequilibrium 



4 




Voltage (mV) [ 

-150 -100 -50 50 100 150 

Voltage (mV) 



FIG. 2: (a) Differential resistance normalized to the normal 
state (taken as a resistance at V ~ 30 mV) for 3 different 
contacts with J?o=45, 43, and 111 fl (curves 1-3, respectively). 
Main gap dV/dl— minima are at about 2.1, 2.6, and 4.86 mV, 
correspondingly. The larger gap for curve 1 equals 7.0 meV. 
Curves 2,3 are shifted vertically for clarity, 
(b) Raw data of the second harmonic signal V2 with modula- 
tion voltages V\ =3.31, 2.78, and 2.5 mV, respectively. The 
numbers of curves are the same as in (a). T = 4.2 K, H = 0. 
Curve 1 are shifted vertically for clarity. Their center of sym- 
metry has (0,0) coordinates. 



FIG. 3: (a) Differential resistances for MgB2-Cu point con- 
tacts, whose normal state spectra are presented in panel b. 
A (1) = 2.55 meV, A (2) = 3.69 meV. T=4.2 K, H=0. (b) 
Second harmonic signals in the normal state. Vi = 2.2 mV, 
T = 41 K. Note the different Y-scales for curves 1 and 2, 
which means that the generation of phonons in the latter 
case is more intensive. For spectrum No.l the vertical ar- 
rows mark the low frequency bumps at about 20 meV. The 
dotted tangent straight lines point to the phonon features on 
curve No. 2. The horizontal bar with double arrows stands for 
thermal smearing, determining the spectral resolution SV. 



phenomena in the superconducting state. These phenom- 
ena are mostly due to destruction of superconductivity in 
the vicinity of the contact, and lead to depression of the 
excess current [Q. It has already been noted, that the 
positions of these nonlinearities are strongly influenced 
by the field and temperature, and should therefore be 
disregarded. 

As to the spectacular electron-phonon spectrum, pre- 
sented in our preliminary report (see Fig. 3, 4 in Ref.Q), 
it is presumably a rare example of a direct contribution 
from the inelastic spectrum of the 2D-band. The mag- 
netic field of 4 T suppresses almost completely the in- 
duced phonon features from the 3D-band (with exception 
of low energy peak). Thus, the proposed tentative con- 



sideration would bring our previous observations in line 
with to what is consistently described here. 

The phonon spectra of a contact with a small value 
of energy gaps are characterized by the presence of low 
frequency phonon peaks. The small peak at energy of 
about 25 meV (Fig. 4, curves 1,2) is visible, where a tiny 
knee exists on the PDOS ||| (see PDOS in Fig. lb). In 
the normal state (at T > T c ), these low frequency peaks 
transform into the S-shape structure in d 2 V '/ 'dl 2 (V)- 
spectra (curve 6), corresponding to the wide minimum 
of dV/dI(V) near zero bias. This low frequency struc- 
ture is hardly due to the remnants of superconducting 
quasigap at T > T c , since it is absent in junction No. 2, 
whose characteristic is shown in Fig. 3. Rather, it could 
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be thought of as strong interaction in the 3D-band with 
very low-frequency excitations, whose origin is not clear 
yet. On increasing the temperature, the low frequency 
peak broadens further like common spectral features do 
in the normal state. 

Since most of the current through the point contact is 
carried by the 3D band with weak EPI, the phonon singu- 
larity in the superconducting state is induced mainly by 
interband scattering from the 2D band with strong EPI. 
The degree of this induction depends on the interband 
coupling. As was noted above, in the clean limit the in- 
terband coupling in MgB2 is low |9 . That was the reason 
for applying the McMillan-proximity-effect model in or- 
der to fit the gap structure in dI/dV(V) characteristic in 
Ref. [Q. This model leads to the small-intensity phonon 
structure. The elastic scattering, which is certainly high 
in thin film and may be even stronger by point contact 
fabrication process, enhances the interband, as well as 
the intraband coupling. Following Arnold ^7j, and as- 
suming the extension of the proximity effect theory in co- 
ordinate space to the two-band model (i.e. "proximity" 
effect in fc-space), we use the corresponding expression 
for the interband coupling, induced by elastic scattering 
in SN sandwiches. In our case, by analogy, the elastic 
scattering leads to the expression 
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ZfA p 1 +i(h/2T)A 2 (E)N 2 (E) 
Zf(E)+i(h/2T)N 2 (E) : 



(3) 



which, like the McMillan expression for h/2r = (not 
shown), derives to Ai(E) ~ A 2 (E), in other words, to 
the large induced phonon structure in the 3D-band. Here 
Zf is the renormalization function in the 3D-band, N 2 (E) 
is the electron DOS in the 2D-band, r is the elastic scat- 
tering time, and Tjv is the interband scattering rate in the 
McMillan model. Thus, as the interband scattering in- 
creases, the intensity of phonon singularity moves to the 
higher energy, characteristic of the 2D-band. This "in- 
ternal proximity effect" explains also the small variation 
in the position (of the order of superconducting energy 
gap) of phonon peaks with various elastic scattering rates 
produced by contact fabrication. 

The influence of magnetic field on the point-contact 
spectra is shown in Fig. 4 (compare curves 1 and 2). It 
is evident that the field can smear out the intensity of 
high energy peaks, which are induced in the 3D-band 
by EPI of the 2D-band. The disappearance of phonon 
peaks at rising field and temperature proves that they 
do not belong to the inelastic backscattering processes, 
which should have the same intensity both in the super- 
conducting and in normal states. It seems more plausible 
that the high energy phonon peaks are due to the elas- 
tic contribution in the excess current (|J) induced by EPI 
from the 2D band, as was already stated above. 

Unfortunately, we cannot estimate the strength of EPI 
by inelastic PC spectroscopy p2||, since we are not able 



FIG. 4: Second harmonic dependences on field and temper- 
ature. The temperature and magnetic field are equal to: 4.2 
K, T; 4.2 K, 4 T; 10 K, 4 T; 20 K, 4 T; 30 K, 4 T; 40 K, T, 
for curves 1-6, respectively. The modulation voltage V\ is 2.2 
mV. A=2.7 meV. The vertical dotted lines mark the phonon 
singularity on the voltage axis. Curves are shifted vertically 
for clarity, their center of symmetry corresponding to (0,0) 
coordinates. 



to destroy the superconductivity completely by magnetic 
field at low enough temperature. The latter is necessary 
to get good resolution. For thin films this experiment 
will hardly promise encouraging results, since the elastic 
mean free path is strongly reduced, and the intensity of 
inelastic backscattering processes is therefore diminished 
by the factor l/d. Note, that for the superconducting 
self-energy effects the short elastic mean free path does 
not suppress essentially the phonon structure j24|. On 
the contrary, for the two-band model elastic scattering 
even enhances the phonon structure due to the " intrinsic 
proximity effect" described above. 
Conclusions 

We summarize our observation as follows. The repro- 
ducible phonon peaks are seen in the superconducting 
state on d 2 V/dI 2 (V), and disappear in the normal state, 
which proves that they are due to the energy dependence 
of superconducting order parameter. The superconduc- 
tivity is presumably due to the 2D-band EPI, which in- 
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duces the A(E) structure in the 3D band. At a small 
value of the superconducting gap the phonon structure 
is weak, and its intensity begins to increase with growing 
A. The robustness against the magnetic field also grows 
notably with increasing A. 

In the normal state the intensity of the inelastic spec- 
trum also correlates with the value of superconducting 
gap. The tendency is observed that the high energy 
phonon peaks become more prominent for a larger gap. 
For a smaller gap, the low-frequency bump appears in the 
spectra. It might be due to EPI with some unknown low- 
frequency excitation Q , but the further work is needed 
to be sure that it is not caused by phonon peaks of the 
nor mal-metal counter electrode (Cu, Ag) 



T |ie final goal of PC spectroscopy is to use sufficiently 
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